Kinetics of chiral phase transition in hot and dense quark matter 
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The kinetics of chiral transitions in hot and dense quark matter is studied via a microscopic frame- 
work (Nambu-Jona-Lasinio model) and a phenomenological model (Ginzburg-Landau free energy). 
We focus on the far-from-equilibrium ordering dynamics subsequent to a quench from the massless 
quark phase to the massive quark phase. The morphology of the ordering system is characterized 
by the scaling of the order-parameter correlation function. The domain growth process obeys the 
Allen-Cahn growth law, L{t) ^ t^^^. We also study the growth of bubble of the stable massive phase 
in a background of the met ast able massive phase. 

PACS numbers: 12.38.Mh, 24.85.+p 



The kinetics of phase transitions, and the ordering pro- 
cess that occurs after a rapid quench in system param- 
eters hke temperature and pressure, plays an important 
and interesting role in different areas of physics [1]. Dur- 
ing the transition, the system develops a spatial struc- 
ture of randomly-distributed domains which evolve with 
time. This ordering process has been extensively studied 
in many condensed matter systems like ferromagnets, bi- 
nary fluids, liquid crystals, etc. In this letter, we study 
an equally fascinating application in high-energy physics, 
i.e., the quark hadron phase transition. 

Heavy ion collision experiments at high energies pro- 
duce hot and dense strongly-interacting matter, and pro- 
vide the opportunity to explore the phase diagram of 
QCD in the plane of temperature (T) and baryon chemi- 
cal potential (/i). Many model studies [3] as well as recent 
lattice studies [3] indicate that at sufficiently large bary- 
onic densities, there is a line of first-order transitions in 
the (/i, T)-plane between a chirally-symmetric phase and 
a broken-symmetric phase. As one moves along the phase 
boundary towards higher T and smaller /i, the first-order 
transition becomes weaker- ending in a second-order crit- 
ical point in the limit of vanishing current quark mass 
or a rapid crossover for nonzero current quark mass 
While the high-T and small-// region of the QCD phase 
diagram has been explored in recent experiments, future 
heavy-ion collision experiments plan to explore the high 
baryon density regime, particularly the region around the 
critical point [5]. 

In this context, the critical dynamics of the chiral tran- 
sition, and its static and dynamical universality proper- 
ties, have attracted much attention recently. Signatures 



of the critical point which had been proposed using on 
equilibrium thermodynamics [6], have been recently in- 
terpreted within a quasi-stationary framework [?] . In the 
present letter, we study the far-from-equilibrium kinetics 
of chiral phase transitions subsequent to a quench from 
a disordered phase (with vanishing quark condensate) to 
the ordered phase. 

To model chiral symmetry breaking in QCD, we use 
the two-flavor Nambu-Jona-Lasinio (NJL) model [8] with 
the Hamiltonian: 

-G - (^7't^)2] . (1) 

To describe the ground state, we take an ansatz with 
quark-antiquark condensates [9] : 

\vac) = exp Jdk qfik)^a ■ k)hi{k)q°/{-k) - h.c. |0). 

(2) 

Here, ^ q are two-component quark and antiquark 
creation operators, and |0) is the perturbative chiral 
vacuum. Further, hi(k) is a variational function re- 
lated to the quark-antiquark condensate as (iptp) = 
-6/(27r)^ Ei=i / sin[2/ii(k)]. This flavor-dependent 
function can be determined by minimising the energy 
at T = or the thermodynamic potential at nonzero 
T and density. Without going into the details which are 
reported elsewhere [10], we write down the expression for 
the thermodynamic potential as 

n{M,p,ti) = -^^^Y, Jdk{ln[l+exp{-/3oJi)]} 
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Here, we have taken vanishing current quark mass, and 
introduce ^ = G[l + l/{4.Nc)]] M -2gps with ps 
(t/St/;) being the scalar density. Further, py = {tp^tp) is the 
vector density, and uj± = (k^ + M^)^/^ ± u where u = 
li—Gpv/Nc. The resulting phase diagram is shown in Fig. 
[TJ We have taken here a three-momentum ultraviolet 
cutoff A = 631 MeV, and the four-fermion coupling G = 
5.074 X 10-^ MeV-^ [12]. With these parameters, the 
vacuum mass of quarks is M :^ 321 MeV. At T = the 
first-order transition takes place at /i 335 MeV. For 
/i = 0, a second-order transition takes place at T c:^ 195 
MeV. The first-order line meets the second-order line at 
the tri-critical point {Ttcp, Ptcp) ^ (74.9, 285.0) MeV. The 
first-order transition is characterized by the existence of 
metastable phases. The limit of metast ability is denoted 
by the dashed lines in Fig. [TJ referred to as the spinodal 
lines. 
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FIG. 1: Phase diagram of the NJL model in the (/i, T)-plane 
for zero current quark mass. A line of first-order transitions 
(green online) meets a line of second-order transitions (blue 
online) at the tricritical point (tcp). We have {iitcp^Ttcp) — 
(285.0, 74.9) MeV. The dashed lines denote the spinodals for 
the first-order transitions. 



Close to the phase boundary, the thermodynamic po- 
tential may be expanded in a power series of the order 
parameter M: 

n{M) = n (o)+^m2+^m^+^m^+- • • = f (m) . (4) 

One can obtain the coefficients in terms of /3 and p by 
Taylor-expanding the thermodynamic potential in Eq. 
([3]). However, we will treat them as phenomenological 
constants which are adjusted to recover the phase dia- 
gram in Fig. [TJ After all, we are interested in studying 
the kinetics of the chiral transition, usually at parameter 
values far from the critical curve in Fig. [1]- where we 
cannot Taylor-expand 11(M, /3,/i). The extrema of the 
potential in Eq. (j4]) are determined by the gap equa- 
tion F\M) = aM ^ bM^ + dM^ = 0. The solutions 
are Mq = 0, and = {-b ± Vb^ - 4.ad)/2d. For 
6 > 0, the transition is second-order, analogous to an M^- 
potential- the stationary points are M = (for a > 0) 



or M = 0, ±M+ (for a < 0). For a < 0, the preferred 
equilibrium state is the one with massive quarks. For 
6 < 0, the solutions of the gap equation are as follows: 
(i) M = for a > (ii) M = 0, ±M+, ±M_ for 

> a > 0, and (iii) M = 0, ±M+ for a < 0. A 
first-order transition takes place at ac = 36^/16 with the 
order parameter jumping discontinuously from M = 
to M = ±M+ = ±(3|6|/4d)^/^ The tricritical point is 
located at btcp = 0, atcp = 0. 

Let us consider the dynamical environment of a heavy- 
ion collision. As long as the evolution is slow compared 
to the typical re-equilibration time, the order parameter 
field will be in local equilibrium. We consider a system 
which is rendered thermodynamically unstable by a rapid 
quench from the disordered phase to the ordered phase. 
The unstable disordered state evolves via the emergence 
and growth of domains rich in the preferred phase [l|, 
13]. The coarsening system is inhomogeneous, and we 
include a surface tension term in the Ginzburg-Landau 
free energy as follows: 



n{M) 



dr 
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The evolution of the system is described by the time- 
dependent Ginzburg-Landau (TDGL) equation 



d_ 
dt 



M{r,t) 



modeling the over-damped relaxational dynamics of 
M(r,t) to the minimum of Vt[M] [l^. Here, F is the 
inverse damping coefficient, and ^(r, t) is the noise term 
satisfying the fluctuation-dissipation relations {0 (r, t)) = 
and {e{Y' ,t')e{Y" ,t")) = 2TT5{y' - Y")5(t' -t"). We 
use the natural scales of order parameter, space and 
time to introduce the dimensionless variables, M = 
./WWl M', r = ^fKM r', t = t'/{r\a\), = 
(r|ap/^/|6|^/^) 0' . Dropping the primes, we have the 
dimensionless TDGL equation: 



-sgn (a) M - sgn (b) 
-\M^^V^M^e(Y,t), (7) 



where A = \a\d/9. For T = 15 MeV, as p takes values 
320 MeV, 330 MeV, 338 MeV and 350 MeV, the corre- 
sponding values of A are 0.065, 0.14, 0.23 and 0.34, which 
we have obtained by fitting the effective potential Eq. (|3|) 
with Eq. dH). 

First, we study the ordering dynamics for 6 > 0. 
We solve the Eq. ([7]) numerically using an Euler- 
discretization scheme. This is implemented on a 3-d lat- 
tice of size N"^ {N = 256), with periodic boundary con- 
ditions in all directions. The dimensionless mesh sizes 



3 



are Ax = 1.0 and At = 0.1, which satisfy the numerical 
stabihty condition. We have further confirmed that the 
spatial mesh size is sufficiently small to resolve the inter- 
face region. The details of the numerics will be reported 
elsewhere llOl . 

t = 20 






FIG. 2: Interface evolution after temperature quench through 
second-order line in Fig. [1] The 3-d snapshots on the left show 
the interfaces (M = 0) at t = 20, 100. They were obtained 
by numerically solving Eq. (0 with a < 0, 6 > 0, A = 0.14. 
The noise amplitude was e = 0.008. The frames on the right 
show a cross-section of the snapshots at 2; = N/2. 




64 




1/2 












16 






(b) 




100 


1000 

t 





FIG. 3: (a) Shows the scaling of correlation function for A = 
0.14 at four different time steps. OJK function (as for usual 
M^-free energy) has good agreement with simulation data, 
(b) The domain size L(t) vs t for A = 0.14. Domain growth 
data is consistent with Allen-Cahn growth law L{t) ^ t^^^ 

In Fig. [21 we show the evolution of a disordered ini- 
tial condition for Eq. ([7]) with a < 0, i.e., a temperature 
quench through the second-order line in Fig. [H The 
initial state consisted of small-amplitude thermal fluctu- 
ations about the massless phase M — 0. The system 
rapidly evolves into domains of the massive phase with 
M and M — M+. The interfaces correspond to 

M = 0- their evolution is shown in the snapshots (frames 
on left) of Fig. O The frames on the right shows the 
interface structure in a cross-section of the snapshots. 
The domains have a characteristic length scale which 
grows with time. The growth process is analogous to 
the coarsening dynamics in the TDGL equation with an 
M^-potential, where growth is driven by kinks with the 
equilibrium profile M{z) = tanh(±2;/A/2). The order- 
parameter correlation function C{r^t) shows dynamical 



scaling or dynamical self- similarity^ C{r^t) = f{r/L)- 
the scaling function f{x) has been calculated by Ohta et 
al. (OJK) [11]. In Fig. [3l^a), we demonstrate that C{r,t) 
for the evolution in Fig. [2] shows the scaling property. 
Further, the domain scale obeys the Allen-Cahn (AC) 
growth law, L{t) ~ t^/^ [see Fig. [3fb)]. The interface 
velocity v ~ dL/dt ~ 1/L, where measure the local 
curvature of the interface. This yields the AC growth 
law. 

Next, let us consider the case with b < 0. In this 
case, a first-order chiral transition occurs for a < ar = 
3|6|V(16(i) (or A < A^ = 3/16). For a < 0, the potential 
has a double well structure and the ordering dynamics is 
analogous to M^-theory i.e., the domain growth scenario 
is similar to Fig. [3l We focus our attention on a quench 
from the disordered state (M = 0) to < A < Ac, corre- 
sponding to a quench between the ffist-order line and Si 
in Fig. [H The massless state is now a metastable state of 
the M^-potential. The chiral transition proceeds via the 
nucleation and growth of droplets of the preferred phase 
(M = ±M+). This nucleation results from large fluc- 
tuations in the initial condition or thermal fluctuations 
during the evolution. In Fig. HI we show the nucleation 
and growth process. At early times {t = 400) the system 
is covered with the M = phase with small droplets of 
the preferred phase. These droplets grow in time and co- 
alesce into domains. The subsequent coarsening of these 
domains is analogous to that in Figs. [2]l3l In the late 
stages, there is no memory of the nucleation dynamics 
which characterized growth during the early stages. 
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FIG. 4: Quenching through first-order phase transition line. 
Evolution of regions with M = (black colored) for A = 0.14 
and for times t = 400 and t = 4000. To the right of the 
each evolution pattern, the order parameter profiles along the 
diagonal are also shown. 

We have also studied the growth of single droplets of 
the preferred phase (M = +M+) in the background of 
the metastable phase (M = 0) [see Fig. Efa)]. The 
droplets have a unique growth velocity v(A), which de- 
pends on the degree of undercooling (A). We have ob- 
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FIG. 5: (a) Shows the growth of single droplet of the preferred 
phase (M = +M+) in the background of the metastable phase 
(M = 0) for A = 0.14 for different time steps. The inner circle 
corresponds to a droplet at time, t = 20. (b) Shows the plot 
of droplet velocity v vs. A. Circles refer to numerical data 
while solid line corresponds to predictions from a phase-plane 
analysis. 



tained v{X) by undertaking a phase-plane analysis of the 
travelling- wave solutions of Eq. ([7]). The droplet inter- 
face corresponds to a saddle connection between the fixed 
points +M+ ^ or -M+ 0. The details of this analy- 
sis will be presented elsewhere [10]. In Fig. [5fb), we plot 
numerical results for v{X) vs. A along with our theoreti- 
cal prediction. Before concluding, we briefly discuss the 
inertial kinetic equation, which is the counterpart of the 
overdamped TDGL equation ([6]). The TDGL equation 
contains first-order time-derivatives and is not Lorentz 
invariant. This is not a serious problem as we are con- 
sidering finite-temperature field theory, where there is 
no Lorentz invar iance. Nevertheless, one can also study 
the kinetics of chiral transitions via an inertial evolution 
equation: 

+ M = -sgn (a) M - sgn (6) 

-AM^ + 6>(r,t). (8) 

In principle, Eq. ([8]) can be obtained from a micro- 
scopic field-theoretic description of the nonequilibrium 
dynamics of the scalar field at finite temperatures [l5| . 
We have also studied the ordering dynamics for the in- 
ertial case, and will present details of our simulations 
elsewhere [16|. Here, we mention the main results of our 
study. The ordering dynamics in this case is similar to 
that in the overdamped case, except that nucleation does 
not have a significant effect even during the early stages 
of evolution. The droplets grow very rapidly and merge 
to form a bicontinueous domain structure characteristic 
of late stage domain growth. The domain growth law is 
again the AC law L{t) ~ t^/^. 

To summarize, we have studied the kinetics of the 
chiral phase transition in QCD. In terms of the quark 
degrees of freedom, the equilibrium phase diagram is 
obtained using the NJL model. An equivalent coarse- 
grained description is obtained from an M^-Ginzburg- 
Landau (GL) free energy. We study the chiral kinetics 
via the corresponding TDGL equation, and consider both 



the overdamped and inertial cases. We study the order- 
ing dynamics resulting from a sudden quench in tem- 
perature through the first-order (I) or second-order (II) 
transition lines in Fig. [TJ For quenches through II and 
deep quenches through I, the massless phase is sponta- 
neously unstable and evolves to the ordered phase via 
spinodal decomposition. For shallow quenches through 
I, the massless phase is metastable and the phase tran- 
sition proceeds via nucleation and growth of droplets. 
The subsequent merger of these droplets results in late- 
stage domain growth analogous to that for the unsta- 
ble case. In all cases, the late-stage growth process ex- 
hibits self-similar dynamical scaling, and the growth law 
is 
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